The intrinsic substrate specificity of the human tyrosine kinome.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Substrate displacement of CK1 C-termini regulates kinase specificity.

CK1 kinases participate in many signaling pathways, and their regulation is of meaningful biological consequence. CK1s autophosphorylate their C-terminal noncatalytic tails, and eliminating these tails increases substrate phosphorylation in vitro, suggesting that the autophosphorylated C-termini act as inhibitory pseudosubstrates. To test this prediction, we comprehensively identified the autophosphorylation sites on Schizosaccharomyces pombe Hhp1 and human CK1ε. Phosphoablating mutations increased Hhp1 and CK1ε activity toward substrates. Peptides corresponding to the C-termini interacted with the kinase domains only when phosphorylated, and substrates competitively inhibited binding of the autophosphorylated tails to the substrate binding grooves. Tail autophosphorylation influenced the catalytic efficiency with which CK1s targeted different substrates, and truncating the tail of CK1δ broadened its linear peptide substrate motif, indicating that tails contribute to substrate specificity as well. Considering autophosphorylation of both T220 in the catalytic domain and C-terminal sites, we propose a displacement specificity model to describe how autophosphorylation modulates substrate specificity for the CK1 family.

Illuminating the Dark Cancer Phosphoproteome Through a Machine-Learned Co-Regulation Map of 26,280 Phosphosites.

Mass spectrometry-based phosphoproteomics offers a comprehensive view of protein phosphorylation, but limited knowledge about the regulation and function of most phosphosites restricts our ability to extract meaningful biological insights from phosphoproteomics data. To address this, we combine machine learning and phosphoproteomic data from 1,195 tumor specimens spanning 11 cancer types to construct CoPheeMap, a network mapping the co-regulation of 26,280 phosphosites. Integrating network features from CoPheeMap into a machine learning model, CoPheeKSA, we achieve superior performance in predicting kinase-substrate associations. CoPheeKSA reveals 24,015 associations between 9,399 phosphosites and 104 serine/threonine kinases, including many unannotated phosphosites and under-studied kinases. We validate the accuracy of these predictions using experimentally determined kinase-substrate specificities. By applying CoPheeMap and CoPheeKSA to phosphosites with high computationally predicted functional significance and cancer-associated phosphosites, we demonstrate the effectiveness of these tools in systematically illuminating phosphosites of interest, revealing dysregulated signaling processes in human cancer, and identifying under-studied kinases as putative therapeutic targets.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

NAK-associated protein 1/NAP1 activates TBK1 to ensure accurate mitosis and cytokinesis.

Subcellular location and activation of Tank Binding Kinase 1 (TBK1) govern precise progression through mitosis. Either loss of activated TBK1 or its sequestration from the centrosomes causes errors in mitosis and growth defects. Yet, what regulates its recruitment and activation on the centrosomes is unknown. We identified that NAK-associated protein 1 (NAP1) is essential for mitosis, binding to and activating TBK1, which both localize to centrosomes. Loss of NAP1 causes several mitotic and cytokinetic defects due to inactivation of TBK1. Our quantitative phosphoproteomics identified numerous TBK1 substrates that are not only confined to the centrosomes but are also associated with microtubules. Substrate motifs analysis indicates that TBK1 acts upstream of other essential cell cycle kinases like Aurora and PAK kinases. We also identified NAP1 as a TBK1 substrate phosphorylating NAP1 at S318 to promote its degradation by the ubiquitin proteasomal system. These data uncover an important distinct function for the NAP1-TBK1 complex during cell division.

Epinephrine inhibits PI3Kα via the Hippo kinases.

The phosphoinositide 3-kinase p110α is an essential mediator of insulin signaling and glucose homeostasis. We interrogated the human serine, threonine, and tyrosine kinome to search for novel regulators of p110α and found that the Hippo kinases phosphorylate p110α at T1061, which inhibits its activity. This inhibitory state corresponds to a conformational change of a membrane-binding domain on p110α, which impairs its ability to engage membranes. In human primary hepatocytes, cancer cell lines, and rodent tissues, activation of the Hippo kinases MST1/2 using forskolin or epinephrine is associated with phosphorylation of T1061 and inhibition of p110α, impairment of downstream insulin signaling, and suppression of glycolysis and glycogen synthesis. These changes are abrogated when MST1/2 are genetically deleted or inhibited with small molecules or if the T1061 is mutated to alanine. Our study defines an inhibitory pathway of PI3K signaling and a link between epinephrine and insulin signaling.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

SARS-CoV-2 hijacks p38ß/MAPK11 to promote virus replication.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the coronavirus disease 2019 (COVID-19) pandemic, drastically modifies infected cells to optimize virus replication. One such modification is the activation of the host p38 mitogen-activated protein kinase (MAPK) pathway, which plays a major role in inflammatory cytokine production, a hallmark of severe COVID-19. We previously demonstrated that inhibition of p38/MAPK activity in SARS-CoV-2-infected cells reduced both cytokine production and viral replication. Here, we combined quantitative genetic screening, genomics, proteomics, and phosphoproteomics to better understand mechanisms underlying the dependence of SARS-CoV-2 on the p38 pathway. We found that p38ß is a critical host factor for SARS-CoV-2 replication in multiple relevant cell lines and that it functions at a step after viral mRNA expression. We identified putative host and viral p38ß substrates in the context of SARS-CoV-2 infection and found that most host substrates have intrinsic antiviral activities. Taken together, this study reveals a unique proviral function for p38ß and supports exploring p38ß inhibitor development as a strategy toward creating a new class of COVID-19 therapies. IMPORTANCE SARS-CoV-2 is the causative agent of the COVID-19 pandemic that has claimed millions of lives since its emergence in 2019. SARS-CoV-2 infection of human cells requires the activity of several cellular pathways for successful replication. One such pathway, the p38 MAPK pathway, is required for virus replication and disease pathogenesis. Here, we applied systems biology approaches to understand how MAPK pathways benefit SARS-CoV-2 replication to inform the development of novel COVID-19 drug therapies.

RBFOX2 modulates a metastatic signature of alternative splicing in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDA) is characterized by aggressive local invasion and metastatic spread, leading to high lethality. Although driver gene mutations during PDA progression are conserved, no specific mutation is correlated with the dissemination of metastases1-3. Here we analysed RNA splicing data of a large cohort of primary and metastatic PDA tumours to identify differentially spliced events that correlate with PDA progression. De novo motif analysis of these events detected enrichment of motifs with high similarity to the RBFOX2 motif. Overexpression of RBFOX2 in a patient-derived xenograft (PDX) metastatic PDA cell line drastically reduced the metastatic potential of these cells in vitro and in vivo, whereas depletion of RBFOX2 in primary pancreatic tumour cell lines increased the metastatic potential of these cells. These findings support the role of RBFOX2 as a potent metastatic suppressor in PDA. RNA-sequencing and splicing analysis of RBFOX2 target genes revealed enrichment of genes in the RHO GTPase pathways, suggesting a role of RBFOX2 splicing activity in cytoskeletal organization and focal adhesion formation. Modulation of RBFOX2-regulated splicing events, such as via myosin phosphatase RHO-interacting protein (MPRIP), is associated with PDA metastases, altered cytoskeletal organization and the induction of focal adhesion formation. Our results implicate the splicing-regulatory function of RBFOX2 as a tumour suppressor in PDA and suggest a therapeutic approach for metastatic PDA.

Depletion of creatine phosphagen energetics with a covalent creatine kinase inhibitor.

Creatine kinases (CKs) provide local ATP production in periods of elevated energetic demand, such as during rapid anabolism and growth. Thus, creatine energetics has emerged as a major metabolic liability in many rapidly proliferating cancers. Whether CKs can be targeted therapeutically is unknown because no potent or selective CK inhibitors have been developed. Here we leverage an active site cysteine present in all CK isoforms to develop a selective covalent inhibitor of creatine phosphagen energetics, CKi. Using deep chemoproteomics, we discover that CKi selectively engages the active site cysteine of CKs in cells. A co-crystal structure of CKi with creatine kinase B indicates active site inhibition that prevents bidirectional phosphotransfer. In cells, CKi and its analogs rapidly and selectively deplete creatine phosphate, and drive toxicity selectively in CK-dependent acute myeloid leukemia. Finally, we use CKi to uncover an essential role for CKs in the regulation of proinflammatory cytokine production in macrophages.

An atlas of substrate specificities for the human serine/threonine kinome.

Protein phosphorylation is one of the most widespread post-translational modifications in biology1,2. With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes3,4. For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible3. Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways.

Host protein kinases required for SARS-CoV-2 nucleocapsid phosphorylation and viral replication.

Multiple coronaviruses have emerged independently in the past 20 years that cause lethal human diseases. Although vaccine development targeting these viruses has been accelerated substantially, there remain patients requiring treatment who cannot be vaccinated or who experience breakthrough infections. Understanding the common host factors necessary for the life cycles of coronaviruses may reveal conserved therapeutic targets. Here, we used the known substrate specificities of mammalian protein kinases to deconvolute the sequence of phosphorylation events mediated by three host protein kinase families (SRPK, GSK-3, and CK1) that coordinately phosphorylate a cluster of serine and threonine residues in the viral N protein, which is required for viral replication. We also showed that loss or inhibition of SRPK1/2, which we propose initiates the N protein phosphorylation cascade, compromised the viral replication cycle. Because these phosphorylation sites are highly conserved across coronaviruses, inhibitors of these protein kinases not only may have therapeutic potential against COVID-19 but also may be broadly useful against coronavirus-mediated diseases.

Evolution of host protease interactions among SARS-CoV-2 variants of concern and related coronaviruses.

Previously, we showed that coagulation factors directly cleave SARS-CoV-2 spike and promote viral entry (Kastenhuber et al., 2022). Here, we show that substitutions in the S1/S2 cleavage site observed in SARS-CoV-2 variants of concern (VOCs) exhibit divergent interactions with host proteases, including factor Xa and furin. Nafamostat remains effective to block coagulation factor-mediated cleavage of variant spike sequences. Furthermore, host protease usage has likely been a selection pressure throughout coronavirus evolution, and we observe convergence of distantly related coronaviruses to attain common host protease interactions, including coagulation factors. Interpretation of genomic surveillance of emerging SARS-CoV-2 variants and future zoonotic spillover is supported by functional characterization of recurrent emerging features.

Coagulation factors directly cleave SARS-CoV-2 spike and enhance viral entry.

Coagulopathy is a significant aspect of morbidity in COVID-19 patients. The clotting cascade is propagated by a series of proteases, including factor Xa and thrombin. While certain host proteases, including TMPRSS2 and furin, are known to be important for cleavage activation of SARS-CoV-2 spike to promote viral entry in the respiratory tract, other proteases may also contribute. Using biochemical and cell-based assays, we demonstrate that factor Xa and thrombin can also directly cleave SARS-CoV-2 spike, enhancing infection at the stage of viral entry. Coagulation factors increased SARS-CoV-2 infection in human lung organoids. A drug-repurposing screen identified a subset of protease inhibitors that promiscuously inhibited spike cleavage by both transmembrane serine proteases and coagulation factors. The mechanism of the protease inhibitors nafamostat and camostat may extend beyond inhibition of TMPRSS2 to coagulation-induced spike cleavage. Anticoagulation is critical in the management of COVID-19, and early intervention could provide collateral benefit by suppressing SARS-CoV-2 viral entry. We propose a model of positive feedback whereby infection-induced hypercoagulation exacerbates SARS-CoV-2 infectivity.

Dietary fructose improves intestinal cell survival and nutrient absorption.

Fructose consumption is linked to the rising incidence of obesity and cancer, which are two of the leading causes of morbidity and mortality globally1,2. Dietary fructose metabolism begins at the epithelium of the small intestine, where fructose is transported by glucose transporter type 5 (GLUT5; encoded by SLC2A5) and phosphorylated by ketohexokinase to form fructose 1-phosphate, which accumulates to high levels in the cell3,4. Although this pathway has been implicated in obesity and tumour promotion, the exact mechanism that drives these pathologies in the intestine remains unclear. Here we show that dietary fructose improves the survival of intestinal cells and increases intestinal villus length in several mouse models. The increase in villus length expands the surface area of the gut and increases nutrient absorption and adiposity in mice that are fed a high-fat diet. In hypoxic intestinal cells, fructose 1-phosphate inhibits the M2 isoform of pyruvate kinase to promote cell survival5-7. Genetic ablation of ketohexokinase or stimulation of pyruvate kinase prevents villus elongation and abolishes the nutrient absorption and tumour growth that are induced by feeding mice with high-fructose corn syrup. The ability of fructose to promote cell survival through an allosteric metabolite thus provides additional insights into the excess adiposity generated by a Western diet, and a compelling explanation for the promotion of tumour growth by high-fructose corn syrup.

Coagulation factors directly cleave SARS-CoV-2 spike and enhance viral entry.

Coagulopathy is a significant aspect of morbidity in COVID-19 patients. The clotting cascade is propagated by a series of proteases, including factor Xa and thrombin. While certain host proteases, including TMPRSS2 and furin, are known to be important for cleavage activation of SARS-CoV-2 spike to promote viral entry in the respiratory tract, other proteases may also contribute. Using biochemical and cell-based assays, we demonstrate that factor Xa and thrombin can also directly cleave SARS-CoV-2 spike, enhancing viral entry. A drug-repurposing screen identified a subset of protease inhibitors that promiscuously inhibited spike cleavage by both transmembrane serine proteases as well as coagulation factors. The mechanism of the protease inhibitors nafamostat and camostat may extend beyond inhibition of TMPRSS2 to coagulation-induced spike cleavage. Anticoagulation is critical in the management of COVID-19, and early intervention could provide collateral benefit by suppressing SARS-CoV-2 viral entry. We propose a model of positive feedback whereby infection-induced hypercoagulation exacerbates SARS-CoV-2 infectivity.

Effects of process factors on the performance of electrochemical disinfection for wastewater in a continuous-flow cell reactor.

Although electrochemical disinfection has been shown to be an effective approach to inactivate bacteria in saline water, the effects of process parameters and reactor design for its application in low-salinity water have not been well understood. In this study, factorial experiments were performed to investigate the direct and confounded effects of applied current (5-20 mA), contact time (2.5-20 min), anode surface area (185-370 cm2), and chloride concentration (50-400 mg L-1) on the disinfection efficiency in fresh water and the secondary effluent of municipal wastewater. An electrochemical disinfection reactor cell with an internal volume of 75 cm3 was designed and fabricated. Residence time distribution analysis showed that the internal mixing of the reactor is similar to that of a dispersed plug-flow reactor. All studied process parameters showed significant effect on the kill efficiency, with the applied current and contact time having the most dominant effect. Although the effect of chloride concentration, which is responsible for electrochemical production of free chlorine in water, is statistically significant, it is not as prominent as those reported for high salinity water. A synergistic effect between chloride concentration and anode surface area was identified, leading to high kill efficiency (99.9%, 3 log kill) at low current density (0.0135 mA cm-2). Response surface modeling results suggested that a scaled-up disinfection reactor can be designed using large anode surface area with long contact time for high chloride water (400 mg L-1) or high current density with short contact time for low chloride water (50 mg L-1). The power requirement of a portable system treating 37.85 m3 day-1 (10,000 gpd) of municipal wastewater was estimated to be 1.9 to 8.3 kW to achieve a 3 log kill, depending on the reactor design.

Identification of SARS-CoV-2 inhibitors using lung and colonic organoids.

There is an urgent need to create novel models using human disease-relevant cells to study severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) biology and to facilitate drug screening. Here, as SARS-CoV-2 primarily infects the respiratory tract, we developed a lung organoid model using human pluripotent stem cells (hPSC-LOs). The hPSC-LOs (particularly alveolar type-II-like cells) are permissive to SARS-CoV-2 infection, and showed robust induction of chemokines following SARS-CoV-2 infection, similar to what is seen in patients with COVID-19. Nearly 25% of these patients also have gastrointestinal manifestations, which are associated with worse COVID-19 outcomes1. We therefore also generated complementary hPSC-derived colonic organoids (hPSC-COs) to explore the response of colonic cells to SARS-CoV-2 infection. We found that multiple colonic cell types, especially enterocytes, express ACE2 and are permissive to SARS-CoV-2 infection. Using hPSC-LOs, we performed a high-throughput screen of drugs approved by the FDA (US Food and Drug Administration) and identified entry inhibitors of SARS-CoV-2, including imatinib, mycophenolic acid and quinacrine dihydrochloride. Treatment at physiologically relevant levels of these drugs significantly inhibited SARS-CoV-2 infection of both hPSC-LOs and hPSC-COs. Together, these data demonstrate that hPSC-LOs and hPSC-COs infected by SARS-CoV-2 can serve as disease models to study SARS-CoV-2 infection and provide a valuable resource for drug screening to identify candidate COVID-19 therapeutics.

Selective inhibition of CDK7 reveals high-confidence targets and new models for TFIIH function in transcription.

CDK7 associates with the 10-subunit TFIIH complex and regulates transcription by phosphorylating the C-terminal domain (CTD) of RNA polymerase II (RNAPII). Few additional CDK7 substrates are known. Here, using the covalent inhibitor SY-351 and quantitative phosphoproteomics, we identified CDK7 kinase substrates in human cells. Among hundreds of high-confidence targets, the vast majority are unique to CDK7 (i.e., distinct from other transcription-associated kinases), with a subset that suggest novel cellular functions. Transcription-associated factors were predominant CDK7 substrates, including SF3B1, U2AF2, and other splicing components. Accordingly, widespread and diverse splicing defects, such as alternative exon inclusion and intron retention, were characterized in CDK7-inhibited cells. Combined with biochemical assays, we establish that CDK7 directly activates other transcription-associated kinases CDK9, CDK12, and CDK13, invoking a "master regulator" role in transcription. We further demonstrate that TFIIH restricts CDK7 kinase function to the RNAPII CTD, whereas other substrates (e.g., SPT5 and SF3B1) are phosphorylated by the three-subunit CDK-activating kinase (CAK; CCNH, MAT1, and CDK7). These results suggest new models for CDK7 function in transcription and implicate CAK dissociation from TFIIH as essential for kinase activation. This straightforward regulatory strategy ensures CDK7 activation is spatially and temporally linked to transcription, and may apply toward other transcription-associated kinases.